Numerical studies of a one-dimensional three-spin spin-glass model 

with long-range interactions 
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We study a p-spin spin-glass model to understand if the finite-temperature glass transition found 
in the mean-field regime of p-spin models, and used to model the behavior of structural glasses, 
persists in the nonmean-field regime. By using a three-spin spin-glass model with long-range power- 
law diluted interactions we are able to continuously tune the (effective) space dimension via the 
exponent of the interactions. Monte Carlo simulations of the spin-glass susceptibility and the two- 
point finite-size correlation length show that deep in the nonmean-field regime, the finite-temperature 
transition is lost whereas this is not the case in the mean-field regime, in agreement with the 
prediction of Moore and Drossel [Phys. Rev. Lett. 89, 217202 (2002)] that three-spin models are 
in the same universality class as an Ising spin glass in a magnetic field. However, slightly in the 
nonmean-field region, we find an apparent transition in the three-spin model, in contrast to results 
for the Ising spin glass in a field. This may indicate that even larger sizes are needed to probe the 
asymptotic behavior in this region. 

PACS numbers: 75.50.Lk, 75.40.Mg, 05.50.+q, 64.60.-i 



I. INTRODUCTION 

There has been considerable interest in relating struc- 
tural glasses to spin glasses because spin-glass models 
are more amenable to analytical and numerical calcula- 
tions than models of interacting atoms. This activity was 
started by Kirkpatrick, Thirumalai, and Wolynesi~— who 
observed a close similarity between the theory for the 
dynamics of p-spin models with p > 2 [at the mean-field 
(MF) level] and mode-coupling theorjsi for the dynamics 
of supercooled liquids. At the mean-field level, the p- 
spin model has two transitions (for a review, see Ref. l5|) . 
There is a dynamical transition at a temperature T = T4, 
also found in mode-coupling theory, below which ergod- 
icity breaking occurs but which is not associated with 
any thermodynamic singularities. In addition, there is a 
transition at T c < Td which does have thermodynamic 
singularities and below which replica symmetry breaking 
(RSB) occurs at the "one-step" level. 5 It is this transi- 
tion which is associated with a possible (ideal) thermo- 
dynamic glass transition of structural glasses, where T c 
corresponds to the Kauzmann temperature 7k £ 

The connection between structural glasses and p-spin 
models is less clear beyond the mean-field level. The dy- 
namical transition at 7^ is an artifact of the mean-field 
limit 7 ^ since it arises from an exponentially large number 
of excited states which trap the system for exponentially 
long times, thereby preventing an infinite system reach- 
ing equilibrium. For a finite-dimensional system, how- 
ever, activation over finite free-energy barriers restores 
ergodicity. Thus the only transition which might occur 
in finite-dimensional p-spin models and structural glasses 
is the thermodynamic transition at T c . 

Even this transition is likely to be significantly differ- 
ent in finite dimensions from mean-field predictions, espe- 



cially for odd p. The reason is that odd-p models violate 
spin-inversion symmetry (Si —Si for all i; Si € {±1}) 
so one might expect that the expectation value of the spin 
would be nonzero at all temperatures T. However, the 
spin average (and hence the spin-glass order parameter) 
is actually zero in mean-field models because of their infi- 
nite connectivity, see, for example, Ref. Nevertheless, 
in any finite-dimensional models, the spin-glass order pa- 
rameter would be nonzero at all T and so any transition 
must be of the replica symmetry breaking type. In fact, 
one of us and Drossel argue that the transition in Jo- 
spin models with odd p is in the same universality class 
as an Ising (p = 2) spin glass in a magnetic field ^ 

Because models with even p have spin-inversion sym- 
metry, which does not seem to have an analog in struc- 
tural glasses, it is natural to take p odd in order to rep- 
resent structural glasses. In the present paper, we study 
numerically whether or not a thermodynamic transition 
occurs in a p — 3 spin glass (and hence presumably also 
in a structural glass) for a range of space dimensions. 

Unfortunately, it is difficult to study spin glasses nu- 
merically in high space dimensions d because the number 
of spins N = L d increases rapidly with linear size L and 
typically one can only study N of order of a few thou- 
sand. Therefore, the range of L is too limited to perform 
a finite-size scaling (FSS) analysis. Recently, it has been 
proposed 1 ^— that one can avoid this difficulty by study- 
ing a model in one dimension in which the interactions 
depend on a power a of the distance. 17 Varying a is anal- 
ogous to varying d in a finite-dimensional model. In this 
paper, we consider values of a corresponding to an effec- 
tive space dimension d c ff both in the mean-field (d e ff > 6) 
and nonmean-field (tf e ff < 6) regions. Our main results 
are that we find a transition in the mean-field region, 
and no transition for a well in the nonmean-field region, 



2 



consistent with our result a 13 ' 16 for the Ising spin glass 
in a magnetic field. However, for a value of a in the 
nonmean-field region, but not far from the critical value 
below which mean-field behavior occurs, we find a tran- 
sition, in contrast to our results for the Ising case. We 
shall discuss possible reasons for this discrepancy. 

The paper is structured as follows: in Sec. [TTl we give 
some theoretical background on the connection between 
the transition in the p > 2 model and that in the Ising 
(p = 2) model in a magnetic field. In Sec. IIII1 we de- 
fine the one-dimensional (ID) three-spin model and de- 
scribe the quantities calculated in the simulations. In 
Sec. IIV1 we briefly give some information on the numer- 
ical method and the parameters of the simulations. Our 
results are presented in Sec. [V] and our conclusions are 
summarized in Sec. I VII 

II. THEORETICAL BACKGROUND 

The field theory associated with p-spin models is a cu- 
bic field theory 1-3 with the following Ginzburg-Landau- 
Wilson Hamiltonian: 



n GLW = / d d r I - £ q%{v) + (V fc/3 (r)) 2 

[ a</3 a</3 

-^Tr g 3(r)-f £^(r)l, (1) 

where q a p is the order parameter and a and (3 are replica 
indices which run from 1 to n, with n — > 0. Terms of order 
q^p and higher have been omitted (and are "irrelevant" 
in the nonmean-field regime). At cubic order there are 
two terms and the ratio of their coefficients R = W2/W1 
plays an important role in the properties of these models 
at the mean-field level. When R > 1, mean-field theory 
predicts 18 that there are two transitions; a dynamical 
transition at Td and a second transition to a state with 
one-step replica symmetry breaking at a lower tempera- 
ture T c . When R < 1, the transitions at Td and T c no 
longer occur; instead there is a single transition to a state 
with full RSB (FRSB). 

Outside the mean-field limit, one-step replica symme- 
try breaking, which occurs in mean field for R > 1, is un- 
stable against thermal fluctuations. 19 As noted in Sec. [TJ 
a FRSB transition, which occurs in mean field for R < 1 , 
is in the same universality class as the Ising spin glass 
in a magnetic fields Therefore, these arguments imply 
that the only possible critical point in finite-dimensional 
p-spin models is in the same universality class as an Ising 
model in a magnetic field. 

A p = 3 model in which the ratio R is less than unity^S 
was numerically studied by Parisi, Picco, and Ritort^ 
who found evidence for a transition. When R < 1, the 
effective field in the Ising spin glass in a field mapping is 



smaller than for R > 1, i.e., the correlation length of the 
system can become very large even if there is no transi- 
tion. When the correlation length becomes on the order 
of the system size, this finite-size effect can be mistaken 
for a genuine phase transition^ It is one of the pur- 
poses of this work to check whether this interpretation 
of the work of Moore and Drossel is correct, by studying 
a p-spin long-range model in one dimension where the 
interplay between the correlation length and the system 
size can be more easily investigated. 



III. MODEL AND OBSERVABLES 




FIG. 1: (Color online) One-dimensional three-spin model. 
The lattice (dashed lines) consists of a two-leg ladder with 
an Ising spin Si at the upper end of the ith rung and an Ising 
spin Ti at the lower end of the rung. An interaction couples 
the two spins at one rung with one of the spins at another 
rung. The solid line shows the interaction involving Si,Ti, 
and Sj. 

We consider a two-leg ladder with Ising spins Si and 
Ti (each take values ±1) on each rung, see Fig. [TJ There 
are L rungs so i = I , . . . , L. Between rungs i and j, 
one can form four combinations of three spins, namely, 
SiTiSj, SiTiTj, SiSjTj, and TiSjTj. With a probability 

Pij ~ r ij 2a > where 

ry = (L/7r)sm(7r|i-i|/L) (2) 

is the geometric distance between the spins arranged on 
a ring, each of these triplets of spins is coupled by an 
independent Gaussian random bond with zero mean 
and standard deviation unity. With a er-dependent prob- 
ability 1 — p^ they are all zero. To avoid the probability 
of placing a bond being larger than 1, a short-distance 
cutoff is applied and thus we take 

Pij = 1 - exp(-C/r^) , (3) 

where the constant C is chosen so that the mean coordi- 
nation number, 

L 

takes a fixed value (z = 6 here). The Hamiltonian is 
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therefore given by 



■CS.SJ, ■ J^I,S,T,) , (5) 



where £y = 1 with probability pij given by Eq. ([3]) and 
zero otherwise. 

We now discuss in detail the correspondence between 
the long-range one-dimensional model in which a is var- 
ied and a short-range spin-glass model in which the di- 
mension d is varied. This correspondence applies quite 
generally for spin-glass models. By varying tr, one can 
tune the model in Eq. ([5]) from the infinite-range to the 
short-range universality classesj 14 ' 16 For < a < 1/2, 
the model is infinite range, in the sense that X^I^'lav 
diverges, and for a = 0, it corresponds to the Viana- 
Bray model^ i.e., a spin glass on a random graph. For 
1/2 < a < 2/3, the model describes a mean- field long- 
range spin glass, corresponding — within the analogy with 
short-range systems — to a short-range model with space 
dimension above the upper critical dimension d > d u = 6. 
For 2/3 < a < 1, the model has nonmean- field critical 
behavior with a finite transition temperature T c . For 
a > 1, the transition temperature is zero. We are inter- 
ested in models which are not infinite range and which 
have a finite T c , i.e., 1/2 < a < 1. 

A rough correspondence between a value of a in the 
long-range ID Ising model and the value of a space di- 
mension d e s in a short-range model can be obtained 
by comparing the scaling of the free-energy density, 
T c t;(T,h,d)~ d , of the d-dimensional system to that in 
the ID long-range system, T C £(T, h.a)^ 1 . When the ex- 
ternal field h is zero, £ ~ \/(T — T c ) v , which gives a 
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matching formula, 

4ff ^SR(4ff) = ^lr(o-). 
A second matching formula^ is 



d e fi = 



2 - 77sR.(4ff) 

2(7-1 



(6) 



(7) 



where ^SR^cff) is the critical exponent rj for the short- 
range model, which is zero in the MF regime. This follows 
from the dependence of £; on h at T c , £ ~ ft, _2 /( d + 2 ~ r '), 
and using the fact that for the long-range system^ 

2 - r; LR = 2ct-1 (MF and non-MF regions). (8) 

Equations ^ and ([7]) agree in the mean-field regime (d > 
6,1/2 < a < 2/3), whereil 



and give 



— !— , u S r=- (MF region) (9) 

Z(7 — 1 Z 



(MF region) . 



(10) 



2(7-1. 

The aforementioned equations also agree to first order 
in 6 — d for d < 6 and at the lower critical dimension. 
Equation ([7]) has the following required properties: (i) 
d c ff — > oo corresponds to a —¥ 1/2, (ii) the upper critical 
dimension d u = 6 corresponds to ct u = 2/3, and (iii) the 
lower critical dimension, which is where di — 2+?7sr(c£z) = 
0, corresponds to a\ = 1. 

To probe the existence of a transition, we compute the 
wave-vector-dependent spin-glass susceptibility given by 



Xso (*) = \ £ ((SiSj)- (Si)(Si}) + ((SiTj)- (Si) (Tj)j \ - (Ti) (S,)) + ((T^) - ( Ti ) (T,-)) ' 



ik (i-j) 



(11) 



where (• • • ) denotes a thermal average and [• • ■ ] av an av- 
erage over the disorder. To avoid bias, each thermal 
average is obtained from a separate copy of the spins. 
Therefore, we simulate four copies at each temperature. 
Note that the spin averages (Si) and (Ti) are nonzero 
even though there is no external field because the inter- 
actions involve three spins and so the model does not 
have spin-inversion symmetry, as discussed in Sec. HI 
The correlation length is given by2%~— 

n 1/(2(7-1) 



6 



1 



Xso(0) 



1 



(12) 



2sin(fc m /2) Lx SG (fc m ) 

where fc m = 2n/L is the smallest nonzero wave vector 
compatible with the boundary conditions. According to 



finite-size scaling,— 

^■=X[L 1 ^(T-T e )], (ct>2/3) 



(13a) 



si Af[L 1 / 3 (T — T c )] , (1/2 < a < 2/3), (13b) 



where i^lr is the correlation length exponent, given in 
the MF region by Eq. ©. Note, from Eq. with 
VSR(dcfi) = 0, which is appropriate for the MF regime, 
and Eq. ((9|), the power of L in Eq. (|13bj) can be re- 
expressed in terms of d g according to 

£l/3(2ff-l) = £<W 6 , (14) 

where the factor of 6 occurs because it is the upper crit- 
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ical dimension d u for spin glasses. The analogous result 
for ferromagnets (for which d u — 4) has been verified 
numerically in Ref. [13] ■ From Eq. (|13|) . if there is a 
transition at T = T c , data for £ L /L (^/i l/LR/3 in the 
mean-field region) should cross at T c for different system 
sizes L. 

We also present data for x SG = X SG (^ = 0)> which has 
the finite-size scaling form 

Xbq =L 2 -^C[L 1 '^{T-T c )\, (<r>2/3), (15a) 
XsG = L^ 3 C[L^ 3 (T - T c )\ , (1/2 < a < 2/3) . 

(15b) 

Hence, curves of X S g/^ 2 ~ ?,lr (Xsg/-^ 3 m the mean-field 
regime) should also intersect. For short-range models, 
Eq. (|15a[) is less useful than Eq. (|13a|) in locating T c be- 
cause it involves an unknown exponent r\. However, for 
long-range models, rj is given by Eq. (j8|) exactly even in 
the nonmean-field regime ) 17 ' 28 i 29 and so Eq. (|15al) is just 
as useful as Eq. (|13a|) in this case<^ 

From now on, all exponents will be those of the long- 
range system so the subscript LR will be suppressed. 

If there are no corrections to scaling, the intersection 
temperatures for all pairs of sizes should be equal to T c . 
However, in practice there are corrections to scaling and 
the intersection temperatures vary with L and only tend 
to a constant for L — > oc. Incorporating the leading cor- 
rection to scaling, which is characterized by a universal 
correction to scaling exponent w, the intersection tem- 
perature of data for, e.g., L and 2L, T*(L,2L), varies 
with L as 



TABLE I: Parameters of the simulations for different values 



of <t. Here iV 8a mp is the number of samples, iV s , 



is the total 



number of Monte Carlo sweeps, T m i n and T max are the lowest 
and highest temperatures simulated, and Nt is the number 
of temperatures. The last column shows the parameter A in 
Eq. ([3]) obtained by fixing z — 6 neighbors on average. 
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N 

1 v samp 


N 

1 ' sweep 
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T 
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Nt 




0.55 


64 


10000 


65536 


1.25 


3.40 


16 


0.95527 


0.55 


128 


10500 


131072 


1.25 


3.40 


16 


0.81746 


0.55 


256 


4400 


524288 


1.25 


3.40 


16 


0.72314 


0.55 


512 


3150 


1048576 


1.55 


3.40 


13 


0.65411 


0.55 


1024 


850 


2097152 


1.55 


3.40 


13 


0.60129 


0.75 


32 


5000 


65536 


0.75 


3.25 


11 


2.02742 


0.75 


64 


5000 


131072 


0.75 


3.25 


11 


1.82345 


0.75 


128 


16300 


524288 


0.75 


3.25 


11 


1.71141 


0.75 


256 


8500 


2097152 


0.75 


3.25 


11 


1.64289 


0.75 


512 


5600 


16777216 


1.00 


3.50 


15 


1.59859 


0.75 


1024 


1000 


33554432 


1.00 


3.50 


15 


1.56903 


0.85 


32 


5000 


65536 


0.25 


4.00 


23 


2.65088 


0.85 


64 


5000 


262144 


0.25 


4.00 


23 


2.47900 


0.85 


128 


4750 


4194304 


0.25 


4.00 


23 


2.39485 


0.85 


256 


3800 


16777216 


0.50 


4.00 


21 


2.34867 



helping to overcome the free-energy barriers found in the 
simulation of glassy systems. 

For the simulations to be in equilibrium, the following 
equality must hold (see Refs. [35j and [16|): 



"-4 



T (l-«) 



(17) 



T*(L,2L)=T C 



(16) 



where A is a nonuniversal amplitude, see Appendix [X] 
and Refs. [3lT[33l Equation (fT6|) is expected to be valid 
in the nonmean-field region, 2/3 < a < 1. Approaching 
the critical value of a = 2/3, one expects ui 0. In 
the mean-field region, 1/2 < a < 2/3, the critical expo- 
nents are known, but we expect corrections to Eq. 1161 as 
discussed in Appendix |A"1 



IV. NUMERICAL METHOD AND 
EQUILIBRATION 

To speed up equilibration, we use the parallel temper- 
ing (exchange) Monte Carlo method. 34 In this approach, 
one simulates Nt copies of the spins with the same inter- 
actions, each at a different temperature between a min- 
imum value T m i n and a maximum value T max . In addi- 
tion to the usual single spin flip moves for each copy, we 
perform global moves in which we interchange the tem- 
peratures of two copies at neighboring temperatures with 
a probability which satisfies the detailed balance condi- 
tion. In this way, the temperature of a particular copy 
performs a random walk between T m ; n and T max , thus 



where U is the energy per rung of the ladder, averaged 
over samples, 



(TiSj) 2 + (TiTtf] (18) 



is the link overlap of a given sample, and N b is the number 
of pairs of connected sites in that sample (i.e., the number 
of nonzero values of ey). In the simulations we keep 
doubling the number of sweeps until Eq. (|17p is satisfied 
within error bars. Note that Eq. (|17p refers to an average 
over samples; the relationship between the energy and 
link overlap is not valid for individual samples. 



V. RESULTS 
A. a = 0.85 

Results for the spin-glass susceptibility divided by 
L 2-n = L 2a-i = L o.7 are ghown in Fi g .[2]f or a = 0.85 

and results for the scaled correlation length are shown in 
Fig. EI The x S g data show no intersections (i.e., no sign 
of a transition). The data for £l/L show an intersection 
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FIG. 2: (Color online) Scaled spin-glass susceptibility for 
a — 0.85, in which 2 — r/ — 2a — 1 = 0.7. According to 
Eq. (|15a[) . the data should intersect at the transition. The 
lack of intersections implies that there is no transition for the 
studied temperature range. 



FIG. 3: (Color online) Scaled spin-glass correlation length for 
a — 0.85. According to Eq. (|13a[l . the data should intersect 
at the transition. Although there is an intersection for the 
smallest pair of size, there is no intersection for the largest 
pair, implying the absence of a transition and in agreement 
with the data in Fig. [2] 



for the smallest pair of sizes, L = 32 and 64 but no in- 
tersection for the largest pair of sizes, L — 128 and 256. 
Hence it appears that for a = 0.85, which is well in the 
nonmean-field regime, there is no transition. Of course, 
we cannot completely exclude a transition at a very low 
temperature. 



B. 



0.75 



Our results for a — 0.75 are shown in Figs. [4] and [5] 
For both Xsg/-^ 2- '' ano - £,l/L, we find nonzero intersec- 
tion temperatures T*(L, 2L) which are plotted in Fig. [5] 
The horizontal axis in Fig. |TJ] is 1/L, and, according to 
Eq. (H~6|) . the data would be a straight line if 1/v+uj = 1. 
Our data are consistent with this but we do not have 
good enough data to obtain a precise value for this expo- 
nent. The main point is that, despite strong corrections 
to scaling, the data for both Xsg/ 1 L 2 ~ v and £,l/L indicate 
a transition with T c in the range from 1.1 to 1.2. 

This is rather surprising since it has been argued^ 
that the transition is in the same universality class as 
the Ising spin glass in a magnetic field, and no transition 
has been found for that model with a — 0.75 in work by 
some of usj 13 i 16 However, corrections to scaling are very 
large (see Figs. HHS]), and so it is plausible that system 
sizes considerably larger than L = 1024 are needed to 
see the true thermodynamic behavior of the three-spin 
model when a \ 2/3, in which case there would be no 




FIG. 4: (Color online) Scaled spin-glass susceptibility for a 
0.75 in which 2-7/ = 2cr-l = 0.5. 



inconsistency with the work of Refs. [lj| and 13 1 
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FIG. 5: (Color online) Scaled spin-glass correlation length for FIG. 7: (Color online) Scaled spin-glass susceptibility for a = 
a = 0.75. 0.55 according to Eq. (|15b[) . The data are consistent with a 

transition at T c ~ 2.1, see also Fig. [9] 



1.5 



S3 



a = 0.75 



O %JL 

■ ^SG 



0.01 



0.02 
1/L 



0.03 



0.04 
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a = 0.55 



1.5 



2.5 



FIG. 6: (Color online) Temperatures where data sets for pairs 
L and 2L intersect for a — 0.75. At large L, the data for both 
X SG and £l/L extrapolate to a value in the range 1.1—1.2. 
This implies that there is a transition at this temperature, 
unless the true asymptotic behavior is only seen at even larger 
sizes. 



FIG. 8: (Color online) Scaled correlation length for a — 0.55. 
The data are consistent with a transition at T c ~ 2.1, see also 
Fig. 03 



C. a = 0.55 



Our results for a = 0.55 (mean-field regime) are shown 
in Figs. [7] and [SJ As discussed in Appendix \K\ the inter- 
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section temperatures in the mean-field regime are given 
by Eq. (|A3I) . For a — 0.55, the exponent 5/3 — 2a is equal 
to 0.57. We therefore plot the intersection temperatures 
against 1/L 057 in Fig. [HI The data strongly suggest that 
there is a transition at T c ~ 2.1 . This result is consistent 
with our earlier results for the Ising spin glass in a mag- 
netic field ) 13 ' 16 where we also found a transition in the 
mean-field region. 



o 



O ML 

■ XsG 



() 



a = 0.55 



which we also find a transition in the mean-field regime, 
this result is seen to be consistent with the predictions of 
Ref. M 

For the case studied that is well in the nonmean-field 
regime, a — 0.85, we find no transition, in agreement 
with our work for the Ising spin glass in a magnetic field. 
This implies that there is no ideal glass transition in three 
dimensions since d — 3 is well below the upper critical 
dimension of d = 6 for models with cubic interactions, 
such as that in Eq. ((T|). 

However, for a = 0.75 the results presented here, which 
indicate a finite transition temperature, appear to be 
at odds with our results for the Ising transition in a 
field; 13 ' 16 where we find no transition. We note, how- 
ever, that Leuzzi et al^ argue that there is a transition 
for this case, based on a nonstandard finite-size scaling 
analysis. In the absence of a transition, the system breaks 
up into domains of size I (Imry-Ma length) which can be 
large at low temperatures, depending on the model. A 
possible explanation of our results for a — 0.75 is that 
£(T — > 0) is greater than the largest system size, namely, 
L = 1024, for the three-spin model, although not for the 
Ising model in a field studied in Ref. If this is the 

case, even larger values of L are needed to determine the 
asymptotic behavior of the three-spin model. 



0.05 0.1 0.15 

1/L 057 

FIG. 9: (Color online) Temperatures where data for pairs L 
and 2L intersect for a = 0.55. At large L, the data for both 
X SG and £l/L extrapolate to a value of approximately 2.1, 
implying that there is a transition at this temperature. 



VI. SUMMARY AND CONCLUSION 
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We have studied the existence of phase transitions in 
a three-spin spin-glass model, that is argued to be an 
appropriate model to describe the (possible) ideal glass 
transition in a supercooled liquid. We have studied three 
values of the parameter a: (i) a = 0.55 (mean-field 
regime), (ii) a = 0.75 (nonmean-field region, but close to 
the mean- field boundary at 2/3), and (hi) a = 0.85 (deep 
inside the nonmean-field regime). Moore and Drosseli^ 
argue that any transition in this model is in the same 
universality class as that of the Ising spin glass in a mag- 
netic field. In particular, the two models should have the 
same critical value of sigma where the transition disap- 
pears (corresponding to the lower critical dimension for 
the short-range case). In other words, if one model has 
a transition the other should have one and vice versa. 

For the mean-field case, a — 0.55, we find a finite- 
temperature transition. Comparing with our previous 
wor k ' 16 for the Ising spin glass in a magnetic field, in 



Appendix A: Size dependence of Intersection 
Temperatures 

According to standard finite-size scaling, the spin-glass 
susceptibility varies near the critical point according to 

XsG (t, L) = L a [f(L b t) + L-»g{LH) + ■ ■ ■} +c + Cl t+- ■ ■ , 

" ^ ; (ai) 

where t = T — T c . The L u term is the leading singu- 
lar correction to scaling and cq is the leading analytic 
correction to scaling. 

Nonmean-field regime. In the nonmean-field 
regime, a c < a < 1 with a c = 2/3, we have a = 2 — rj = 
2a — l and b = Ijv. (In this section all exponents refer 
to the long-range universality class.) We use Eq. (|A1[) 
to calculate the temperature T*(L,2L) where data for 
X SG /L a for sizes L and 2L intersect. Expanding f(x) 
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to first order in x, replacing g(x) by 5(0), and assuming 
that a > oj (which is certainly true near a = 2/3, where 
cj — s- 0) we recover Eq. (fl"6|) . 

Mean-field regime. Curiously, the situation in the 
mean-field regime, 1/2 < a < 2/3, is more compli- 
cated. First of all, the exponents a and b are inde- 
pendent of u 27 i 36 ' 37 and take the value at er c for all 
1/2 < a < a c , i.e. a = b = 1/3. Second, although 
the _L 2tT_1 term is replaced as the largest term by an 
L 1 / 3 term (due to the presence of a "dangerous irrele- 
vant variable," cf. Rcfs. [27.36.37]) we expect this term to 
not disappear but rather become a correction to scaling. 
Hence, we replace Eq. (jAll) by 



XsG (t, L) = L 1 / 3 \f(L^t) + L-«g{l}lH) + 



daL^^hg^H) + c + • • • (1/2 < a < 2/3). 



The correction exponent uj can be obtained in the mean- 
field regime from the work of Kotliar et alw^ and is 3 ^ 
given by uj = 2 — 3c. 

For a < a c , we find that the L 2cr_1 term gives the 
leading correction in Eq. (IA2[) and, as a result, Eq. ([To]) 
is replaced by 



T*(L,2L) = T C - 



A' 



£5/3 -2cr ' 



(1/2 < cr < 2/3). (A3) 



(A2) 



To determine the intersection temperatures of the cor- 
relation length we also need the FSS scaling form for 
% SG (A: m ). We find that there is an additional correction 
which dominates for a < 7/12 = 0.5833. However, for 
the value a = 0.55 used in the simulations, the resulting 
difference from Eq. (IA3|) is very small and therefore we 
neglect it. 
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